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Abstract 

We study the rare decays B + D^+K^°, which can occur only via 
annihilation type diagrams in the standard model. We calculate all of the 
four modes, B — > PP, VP, PV, VV, in the framework of perturbative QCD 
approach and give the branching ratios of the order about 10~ 6 . 

I. INTRODUCTION 

More and more data of B decays are being collected at the two B factories, Belle and 
BaBar. The original approach to non-leptonic B decays based on the factorization approach 
(FA) [1], which succeeded in calculating the branching ratios of many decays [2]. FA is 
a simple method, by which non-factorizable and annihilation contributions are neglected. 
Although calculations are easy in FA, it suffers the problems of scale, infrared-cutoff and 
gauge dependence [3]. Especially it is difficult to explain some observed branching ratios 
of B decays, such as B -> J/ipK 1 -^ [4]. To improve the theoretical application [5] and 
understand why the simple FA works so well [6,7], some methods have been brought forward 
and developed. One of them is the perturbative QCD (PQCD) approach developed by 
Brodsky and Lepage [8,9], under which we can calculate the annihilation diagrams as well 
as the factorizable and non-factorizable diagrams. 

It is consistent to calculate branching ratios of B decays in PQCD approach, as we will 
explain its framework in the next section. It has been applied in the non-leptonic B decays 
[6,7,10,11] successfully. In the case of B + -> £>M+.K-M° decays, which is a kind of pure 
annihilation type decays, the physics picture of PQCD is as follows, shown in Fig. LAW 
boson exchange causes bu — > sc, and the dd quarks are produced from a gluon. In the rest 
frame of the B meson, the d and d quarks in the ZjW+if W° mesons each has momentum 
0(Mb/4), so the gluon producing them has q 2 ~ 0(M§/4). It is a hard gluon according to 
the mass of B meson. So we can perturbatively treat these decays and use PQCD approach 
like other pure annihilation type B decays [12] . 
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In the next section, we explain the framework of PQCD briefly. In section 3, we give 
the analytic formulae for the decays B + — > ZjM+lfM . In section 4, we show the numerical 
results and theoretical errors of the four modes respectively. Finally, we draw a conclusion 
in section 5. 

II. FRAMEWORK 

The PQCD approach divides the process into hard components, which are treated by 
perturbative theory, and non-perturbative components, which are put into the hadron wave 
function. The hadron wave function can be extracted from experimental data or calcu- 
lated by QCD sum rules method. The decay amplitude can be conceptually written as the 
convolution 

Amplitude ~ Jd%d%d% Tr[C(t)$ B (ki)$D|(k2)$K(k 3 )H(k 1 ,k 2 ,k3,t)], (1) 

where fcj's are momenta of light quarks included in each meson, and Tr denotes the trace over 
Dirac and color indices. The hard components comprise hard part (H) and harder dynamics 
(C). H(t) describes the four quark operator and the spectator quark connected by a hard 
gluon. It can be perturbatively calculated, since it includes the hard dynamics characterized 
by the scale t, where t ~ O (M B /2) for B + — »■ D^ + K^° decays, and the hard gluon's q 2 is 
of the order t 2 . C(t) is the Wilson coefficient which results from the radiative corrections at 
short distance. In the above convolution, C(t) includes the harder dynamics at a larger scale 
than the Mb scale and describes the evolution of local four-Fermi operators from Mw down 
to the scale t. The wave function <3>m denotes the non-perturbative components, which is 
independent of the specific processes and removes the infrared cut off dependence in PQCD 
approach. 

According to the conservation of four-momentum, we can obtain £>W+ and meson s 
energy and momenta in the rest frame of B meson, 

E 1 = M B , |Pi|=0, 

E 2 = ^{l + r 2 2 -rll £ 3 = ^(l-r 2 + r 2 ), 

|P 2 | = |P 3 | = ^ V 1 " 2 r 2 2 + r 2 4 - 2 r 3 2 - 2 r 2 2 r 3 2 + r 3 4 , (2) 

where the subscript (1,2,3) denote B + , £)(*)+ and meson respectively, and r 2 = 

M D (*)+ /Mb, t 3 = M K (*)o/Mb- It is convenient to assume that the D* + (K*°) meson moves in 
the plus (minus) z direction carrying the momentum P 2 (-£3)- The longitudinal polarization 
vectors of the D* + and K*° are given as 

e 2L = t^-( \J(l ~ r 3 2 ) 2 - 2 r 2 2 (1 + r 3 2 ) + r 2 4 , 0, 0, 1 + r 2 2 - r 3 2 ), 
Zr 2 

esL = ^-( V(l - r 3 2 ) 2 - 2 r 2 2 (1 + r 3 2 ) + r 2 4 , 0, 0, -1 + r 2 2 - r 3 2 ), (3) 
zr 3 
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which satisfy the normalization e 2L = e 3L = — 1 and the orthogonality e 2 L • P2 = ^3L • P3 — 0. 
For simplicity, we use the light-cone coordinate 1 to describe the meson's momenta in the rest 
frame of B meson. After deducing the analytic formulae of amplitudes, we ignore the terms 
proportional to r\ ~ 0.15 or r\ ~ 0.04 2 . Equivalently we ignore the terms proportional to r 2 
in B + , and K^ meson's momenta and longitudinal polarization vectors. Therefore 

eqs.(2,3) in the light-cone coordinate are corresponding to 

Pi = ^(U, Or), P 2 = ^f(l-r 3 V 2 2 ,0 T ), P 3 = ^|(r 2 ,l-r 2 ,0 T ), (4) 

e 2 L = ^(l-r 3 2 ,-r 2 ,0 T ), e 3L = -^(-r 2 , 1 - r 2 , T ), (5) 
V2r 2 V2r 3 

respectively. The transverse polarization vectors can be adapted directly as e 2 r = (0, 0, It), 
£3t = (0, 0, It)- We denote the light (anti-)quark momenta in B + , and if(*)° mesons as 

ki = (xiPi, 0, kir), k 2 = (x 2 P 2 , 0) k 2 r), and k 3 = (0, X3P3 , k 3 T) respectively. Integrating 
eq. (1) over fcf , fc^", and fcj, we obtain 

Amplitude ~ Jdxidx 2 dx 3 bidbib 2 db 2 b 3 db 3 , (6) 

Tr[C(t)$ B (xi, b 1 )$ Ds (x 2 , b 2 )$ K (x 3 , b 3 )H(x ; , b ; , t)S t (x ; ) e~ s % (7) 

where bi is the conjugate space coordinate of kiT, and t is the largest energy scale in H, as a 
function in terms of Xi and bi. The large logarithms ln(My^/t) resulting from QCD radiative 
corrections to four quark operators are absorbed into the Wilson coefficients C(t). The 
inclusion of k? brings in one kind of large logarithms ln 2 (P6) from the overlap of collinear and 
soft gluon corrections, P denoting the dominant light-cone component of meson momentum. 
The other kind of large logarithms ln(tb) derives from the renormalization of the ultraviolet 
divergences. These two kinds of large logarithms are summed and lead to Sudakov form 
factor, e~ s( - f \ It suppresses the soft dynamics effectively [13]. The large double logarithms 
ln 2 a;j are summed by the threshold resummation [14], and they lead to S t (xi) which smears 
the end-point singularities on X{. From the brief analysis above, it can be seen that PQCD 
is a consistent approach. 



III. ANALYTIC FORMULAE 

A. The wave functions 

We use the wave functions $M,a/3 decomposed in terms of spin structure. The coming B 
meson and outgoing £)(*)+, if (*)° are as follows: 



1 We use the light-cone coordinate in the convention, where = ^(p° ± p 3 ) and pr = (p 1 ^ 2 )- 

2 This approximation is also adapted in deriving meson wave functions. So it is consistent to take 
eqs.(4,5). Moreover we ignore the mass of light pseudo-scalar meson K°. 
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$ B (x,b) 



[(A7 5 )+M B75 ] M^b), 



(8) 



$ D (x,b) = 



[ 75 (A + M D )]0 D (x,6), 



(9) 



<$>d* (x,b) 



[ftf>2 + M D .)]</> D .(x,b), 



(10) 



75 ft<t>i{x, b) + M 0Kl5 (j)^(x, b) 
+M 0Kl5 (j6 /i-l)(j) T K (x,b) 



(11) 



$K*L(x,b) = 



M K * fa L <t>K*(x,b)+ izL f3^ K *(x,b) 
+M K *I(f) s K «(x,b) 



(12) 



$k*t(x, b) 



+ 



M K * faT<t>K*{x)+ far fo<t> T K*{x) 



P 3 -n 



(13) 



where N c = 3 is color's degree of freedom, and M k = M\j{rn u + m s ), v = (0, 1, Ot) oc P3, 
n = (1,0, Ot), e 0123 = 1- The subscripts L and T denote the wave functions corresponding 
to the longitudinally and transversely polarized K* mesons. 



B. The effective Hamiltonian 



The effective Hamiltonian for decay B + — > ipW+ifM at a scale lower than Mw is given 
by [15] 



G 



H cS = ^V: b V cs [CMOM + C 2 (ii)0 2 Qi)] , 
Oi = (Is) V -a(cu)v-a, 2 = (bu) V - A (cs) V - A , 



(14) 
(15) 



where Ci^/i) are Wilson coefficients at renormalization scale //, and summation in SU(3) C 
color's index a and chiral projection, Sa?a7^(l — 7s) 9a are abbreviated to (qq')v-A- The 
lowest order diagrams of B° — > D*~K + are drawn in Fig. 1. We will choose |V^ S | = 
0.996 ± 0.013, |Kb| = (3.6 ± 0.7) x 10~ 3 [16]. There is no CP violation in the decays, since 
only one kind of CKM phase appears in the processes. Therefore the decay width for CP 
conjugated mode, B~ — > D^~K^°, equals to B + — > (*)° respectively. 
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C. The decay width 



The total decay amplitude for each mode or helicity state of B + — > D^ + K W° is written 

as 

A = f B F + M, (16) 

where f B is the decay constant of B meson, and the overall factor is included in the decay 
width with the kinematics factor. F(M) stands for the amplitude of (non-) factorizable 
annihilation diagrams in Fig. la,b (c,d). We exhibit their explicit expressions and subscripts 
of F and M according to the modes and helicity states respectively in the appendix. The 
decay width for each mode of these decays is given as 

r = ^ (1_ * |r - U|1 ' (17) 

where the subscript a denotes the helicity states of the two vector mesons with L(T(1,2)) 
standing for the longitudinal (transverse) component in the case of B + — > D* + K*° decay, 
as shown in the appendix. 



IV. NUMERICAL RESULTS 



In the numerical analysis, we adopt the B meson wave function as [6,7] 

Ml x 2 1 



4> B (x, b) = N B x (1 — x) exp 



2uj 



(18) 



with the shape parameter Ub and the normalization constant N B being related to the decay 
constant f B by normalization 



o 1 dx M x,b = 0) = J^=, 



(19) 



which is also right for meson, i.e. M = B, D^*\ 

For meson wave function, we use two types. The first kind [17] is 



(-0 (f>D(*)( x ) 



■f D (*)x(l -x){l + a D( *)(l - 2x)}exp 



(20) 



in which the last term, exp — |(cj D (,)6) 2 , derived from the distribution. By taking same 
parameters, we neglect the difference between the D and D* mesons wave functions, since 
the c quark is much heavier than the d quark, and the mass difference between two mesons 
is little. The second kind [18] is 



(//) <p DM (x) 



■f D (*)x(l - x){l + a D (*)(l - 2x)}, 



(21) 
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which is fitted from the measured B — > D^lv decay spectrum at large recoil. The absence 
of the last term like (J) is due to the insufficiency of the experiment data. 

The K^> wave functions [19,20] we adopt are calculated by QCD sum rules. To abridge 
the context, we list them and the corresponding parameters in the appendix. 

The other input parameters are listed below: 

f B = 190 MeV , f D = 240 MeV , f D * = 240 MeV , 

f K = 160 MeV , f K . = 200 MeV , /£. = 160 MeV , w flW = 0.2 GeV, (22) 
M K = 1-60 GeV , u b = 0.4 GeV , a DW = 0.3 , C D = 0.8 , C D * = 0.7 , (23) 
M B = 5.279 GeV , M b = 4.8 GeV , M D = 1.869 GeV , M D * = 2.010 GeV , 
M t = 170 GeV , M w = 80.4 GeV , t b ± = 1.674 x 10~ 12 s , G F = 1.16639 x 10~ 5 GeV" 2 , (24) 

where the Fermi coupling constant Gp, the masses and life times of particles refer to [21]. 

With the analytic formulae and parameters above, we get the branching ratios of B + — > 
(*)° shown in Table I, II, III, IV for two kinds of wave functions respectively. The 
magnitude according to wave function I is about 60 percent of the one corresponding 
to wave function II. The difference can tell the correct D^> wave function by the 
experiment data in the future. 

For each mode of B + — > £)(*)+ decays, the contribution of the factorizable and 
non-factorizable annihilation diagrams is the same order, although F is proportional to 
Wilson coefficient Ci + Ci/3, which is O (1), and non-factorizable annihilation diagram 
contribution is proportional to Ci, which is about 30 percent of C 2 + Ci/3. Since the 
counteraction influence between Fig. Va,b of F is heavier than that between Fig. Vc,d of 
M by the reason of the more similar propagators in Fig. Va,b. The magnitude comparison 
can bee seen directly from Table I, II, III, IV. 
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TABLES 





B -► DK 


B -► d*k 


5 -► DK* 


f B F 


-1.56 + 1.21* 


-1.86 + 2.20* 


1.17- 1.65* 


M 


1.03 + 1.29* 


-1.01-0.78* 


-1.75-0.89* 


A 


-0.52 + 2.50*' 


-2.87+1.41* 


-0.58 - 2.54*' 


Br(l(T 6 ) 


0.93 


1.42 


0.96 



TABLE I. The branching ratios of the three decay modes and amplitudes (10~ 2 Gey) in terms 
of the factorizable, non-factorizable diagrams and the sum of them according to the £)W wave 
function I. 





(B - D*K*) L 


(B - L>*K*) T1 


(5 - L>*K*) T2 


f B F 
M 


-0.17-3.44* 
1.73 + 0.35* 


1.58-3.47* 
-0.09 - 0.41* 


-0.37 - 0.22*' 
0.009 + 0.005*' 


A 


1.56 - 3.08*' 


1.49-3.87* 


-0.36 - 0.21*' 


Br(10- 


" 6 ) 


1.67 


2.40 


0.02 


Total Br(10" 6 ) 


4.09 


TABLE II. The branching ratios of B — ► D*K* decay and helicity amplitudes (W~ 2 GeV) in 
terms of the factorizable, non-factorizable diagrams and the sum of them according to the D* wave 
function I. 






B -► L>K 


5 -► D*K 


B -► DK* 


JbF 
M 




-2.38 + 1.56* 
1.37 + 1.49*' 


-2.47 + 2.94*' 
-1.24- 0.87*' 


1.50 - 2.32*' 
-2.18 - 1.01*' 


A 




-1.01 + 3.05*' 


-3.71 + 2.07*' 


-0.68 - 3.32*' 


Br(10~ 6 ) 




1.47 


2.52 


1.64 



TABLE III. The branching ratios of the four decay modes and amplitudes (W~ 2 GeV) in terms 
of the factorizable, non-factorizable diagrams and the sum of them according to the wave 
function II. 
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(B -► D*K*) L 


(B -► D*K*) T1 


(5 -► L>*K*) T2 


f B F 


-0.38 - 4.6H 


2.02 - 4.62i 


-0.44 - 0.22i 


M 


2.03 + 0.37i 


-0.14-0.49z 


0.01 + 0.006i 


A 


1.65 - 4.24i 


1.88-5.11* 


-0.43 - 0.22i 


Br(l(T 6 ) 


2.89 


4.12 


0.03 


Total Br(10" 6 ) 


7.04 



TABLE IV. The branching ratios of B — > D*K* decay and helicity amplitudes (10 _2 Gey) in 
terms of the factorizable, non-factorizable diagrams and the sum of them according to the D* wave 
function II. 
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From Table II, IV, we can see \Ati\ > |Al| ^> |Ar 2 | in the case of B — > VV mode. 
There are two questions worthy of asking, why is |Ar 2 | so little? why are \At\ \ and \Ai\ the 
same order, though \A T i\ is suppressed at least by the term r 2 (r = r 2 or r 3 )? According 
to the amplitudes of F T2 and M T2 , the contribution of the twist 2 wave function § T K « is 
absent, and the coefficients corresponding to the twist 3 wave functions (j) v K , and (j) a K * are 
just opposite and counteract each other heavily. Therefore the value of At2 is too little to 
consider. To answer the second question, we should note that r 2 is not a serious suppression 
term, especially when r 2 times 2, 2r 2 ~ 1, like the term in F?i and M^. In the case of F?i 
and M T1 , all of the signs of the sub- amplitudes corresponding to the two twist 3 K* wave 
functions are same, and the terms in the front of the twist 2 wave function 0^* do not suffer 
the heavy suppression of r^. On the other hand, in Fl (Ml) the seemly main contribution 
of the twist 2 wave function <pK* is offset by the opposite coefficients in Fig. Va,b (c,d). 
Moreover in Fig. Va the signs of the coefficients corresponding to the twist 3 wave function 
<^ K * and (p s K , are different. For the reasons above, \A T1 \ and \A L \ are the same order. 

It should be stressed that there is no arbitrary parameter in our calculation, but we only 
know the magnitude of each up to a range. In Table V,VI we show the sensitivity of the 
branching ratios to 30% change of the parameters in eq. (23) according to the two kinds of 
Z)W wave functions respectively. Since the M 0K and oV s uncertainty influences the results 
very much, we will limit them to a more appropriate extent. According to [19], 

1.4 GeV < M 0K < 1.8 GeV, (25) 

the branching ratios are 

Br(B + D+ K», = { ; : ^| * i u m 

Br( B+ d . + k») = { x io- i u (27) 

where /(//) stands for the result for /(//) kind of wave function. From the B — > K 
transition form factor f+(0), we can limit the appropriate extent of u>b- f+(0) calculated 
from PQCD at m K = 1-6 GeV is consistent with f+(0) by QCD sum rules [19], when 

0.35 GeV < u b < 0.46 GeV. (28) 

In the above range, the branching ratios are 

Br(B+ D + K») = { J^S * \ _ 6 (29) 



Br(B+ -> D*+K°) = 



1.42; 

2.52; 



-o.oo 
-o.oo 
-o.oo 
-o.oi 



10 
10 



(30) 



BrfB" 



D + K 



0.96" 
1-6418 



o.io 

0.09 x 
16 



17 



X 



10~ 6 

io- 6 



I 

II 



(31) 
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twr+ _^ n*+x<'*o\ _ J 4-09to;[!g x 10 6 / 

K ) — | 7 04 +o.07 x 10 -6 JJ 



lvl OK 


RrfR — v FlK"! 


RrfR — v D*K"l 


RrfR — v DK*"l 


RrfR — v 

JDI ID * I_I IV J 


1.1Z 


U. <y 


n ei 
U.ol 






1 Kn 


U.yo 


1 /( 
1.4Z 






2.08 


1.09 


2.22 








Br(B -► DK) 


Br(B -»• D*K) 


Br(B -► DK*) 


Br(B -► D*K*) 


0.28 


1.16 


1.41 


1.21 


4.25 


0.40 


0.93 


1.42 


0.96 


4.09 


0.52 


0.75 


1.42 


0.79 


3.99 


a D M 


Br(B -► DK) 


Br(B -»• D*K) 


Br(B -► DK*) 


Br(B -► D*K*) 


0.21 


0.88 


1.32 


0.90 


3.79 


0.30 


0.93 


1.42 


0.96 


4.09 


0.39 


0.99 


1.54 


1.03 


4.40 



TABLE V. The sensitivity of the branching ratio (10 6 ) to 30% extent of parameters in terms 
of the four modes of the B + — ► D^ + K^° decays according to the wave function I. 



M 0K 


Br(B -► DK) 


Br(B -»• D*K) 


Br(B -► DK*) 


Br(B -► D*K*) 


1.12 


1.17 


1.41 






1.60 


1.47 


2.52 






2.08 


1.80 


3.95 








Br(B -► DK) 


Br(B -► D*K) 


Br(B -► DK*) 


Br(B -► D*K*) 


0.28 


1.86 


2.48 


2.06 


7.22 


0.40 


1.47 


2.52 


1.64 


7.04 


0.52 


1.20 


2.51 


1.34 


6.92 


Cd 


Br(B -► DK) 


Br(B -► D*K) 


Br(B -► DK*) 


Br(B -► D*K*) 


0.56 


1.26 




1.38 




0.80 


1.47 




1.64 




1.04 


1.70 




1.92 




c D * 


Br(B -► DK) 


Br(B -► D*K) 


Br(B -► DK*) 


Br(B -► D*K*) 


0.49 




2.13 




5.99 


0.70 




2.52 




7.04 


0.91 




2.95 




7.50 



TABLE VI. The sensitivity of the branching ratio (10 6 ) to 30% extent of parameters in terms 
of the four modes of the B + — > Z)(*)+Jf (*)° decays according to the wave function II. 
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Besides the Perturbative annihilation contribution above, there is also contribution 
from the final state interaction (FSI) in hadronic level, such as B + — > £)(*) i^(*)+ then 
£)(*)o^(*)+ £)(*)+ Based on the argument of color transparency [9,22], FSI effects 
may not be important in the two-body B decays. So we suppose that the dominant con- 
tribution is what we calculated above. The hypothesis is consistent with the argument in 
[6,23]. 

V. CONCLUSION 

In this paper, we study the four modes of B + — > D^ + K^° decays. Based on the 
consistent PQCD framework, we predict the branching ratios of these pure annihilation 
type decays of the order of 1CT 6 , and show the theoretical errors. Such results can be 
measured in the two B factories in the future. 
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APPENDIX A: APPENDIX 

1. The (non-)factorizable amplitude 

At first order of a s , we get the analytic formulae of the (non-)factorizable amplitude for 
each mode or helicity state listed below. We neglect the small term x\ in the numerators 
of the hard part of M, since the B meson wave function in eq. (18) have a sharp peak at 
the small x region, 0{~h/M B ), where A = M B — M b . It should be noticed that we do not 
employ this approximation to the denominators of the propagator which are sensitive to x 1 . 
Because X\ there behaves as a cut-off. 

a. B+ -► D+K° decay 
The amplitude for the factorizable annihilation diagrams in Fig. la,b is given as 

Fx = 16ttC f M b / dx 2 dx 3 / b 2 db 2 b 3 db 3 (f) D (x 2 ,b 2 ) 
Jo Jo 

x {(x 3 - 2x 3 r 2 2 - r 2 2 )<j)£{x 3 , 6 3 ) + r 2 r K (1 + 2x 3 ) 0^(x 3 , 63) - r 2 r K {l - 2x 3 )4>l{x 3 , b 3 )} (Al) 
E f(tl)h a (x 2 , x 3 , b 2 , b 3 ) + {{r- 2 2 - l)x 2 (f)^(x 3 , 63) - 2r 2 r^(l + x 2 )(j)^(x 3 , b 3 )}E f (t 2 a )h a (x 3 , x 2 , b 3 , b 2 

The amplitude for the non-factorizable annihilation diagrams in Fig. la,b is obtained as 
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Mi = 



j rl roo 

7= _ 647tCfM| dx 1 dx 2 dx 3 b 1 dbib 2 db 2 bi)4> D (x 2 , b 2 ) 

yZJMc Jo Jo 



x 



{(x 3 + (x 2 - 2:r 3 - I)r2)0^(x3, 6 2 ) + r 2 r K (2 + x 2 + x 3 ) 0^(x 3 , b 2 ) - r 2 r K (x 2 - x 3 ) 4> T K (x 3 , b 2 )} 
Em(t 1 m )h£ ) (x 1 ,x 2 , x 3 ,h, b 2 ) 

-{x 2 4>k(x 3 , b 2 ) + r 2 r K (x 2 + x 3 ) 0^(x 3 , b 2 ) + r 2 r K (x 2 - x 3 ) <Pk(x 3 , b 2 )}E m (t 2 m )h (2) (x 1 , x 2 , x 3 , b ± , b 2 ) 

where C F = 4/3 is the group factor of SU(3) C gauge group, and r K = M 0K /M B , and the 
functions Ef, E m , t^' 2 , h a are given in the appendix A. 3. 



b. B + -► D*+K° decay 

rl roo 

F 2 = -1<6tiC f M 2 b / dx 2 dx 3 / b 2 db 2 b 3 db 3 (fi D *(x 2 ,b 2 ) 
Jo Jo 



x 



{(x 3 - 2x 3 r\ + rl)(p^(x 3 , b 3 ) + r 2 r K (f K {x 3 , b 3 ) - r 2 r K (j) T K {x 3 , b 3 )}E f (t 1 a )h a (x 2 , x 3 , b 2 , b 3 ) 



-{(1 - rj)x 2 (p^(x 3 , b 3 ) - 2r 2 r K {\ - x 2 )(p^(x 3 , b 3 )}E f (t 2 a )h a (x 3 , x 2 , b 3 , b 2 ) 



I rl roo 

M 2 = —==QAttC f MI / dxidx 2 dx 3 \ bxdb x b 2 db 2 4>B(xi,b 1 )(j) D *(x 2 ,b 2 ) 
V2iV c Jo Jo 

x \{(x 3 + (l-x 2 - 2x 3 )rl)(p^(x 3 , b 2 ) + r 2 r K (x 3 - x 2 ) 0^(x 3 , b 2 ) 

-r 2 r K (2-x 2 - x 3 ) (f> T K {x 3 , b 2 )}E m (t 1 m )h^\x 1 , x 2 , x 3 , b ± , b 2 ) 
-{(1 - 2rl)x 2 (px(x 3 , b 2 ) + r 2 r K (x 2 - x 3 ) 0^(x 3 , b 2 ) 
+r 2 r K (x 2 + x 3 ) 4>k(x 3 , b 2 )}E m (t 2 m )h {2) (xi, x 2 , x 3 , b ± , b 2 ) . 



(A3) 



(A4) 



c. B+ -► D + K*° decay 

rl roo 

F 3 = -IQtiC f M% I dx 2 dx 3 f b 2 db 2 b 3 db 3 (f)D(x 2 ,b 2 ) 
Jo Jo 

x [{(x 3 - 2x 3 r\ - rl)(j) K *(x 3 ,b 3 ) +r 2 r 3 (1 + 2x 3 ) 4> s K *{x 3 , b 3 ) - r 2 r 3 (l - 2x 3 )(f) t Kt (x 3 ,b 3 )} 
Ef{t l a )h a {x 2 ,x 3 ,b 2 ,b 3 ) 

-{(1 - r 2 2 )x 2 (t)K*{x 3 M) + 2r 2 r 3 (l + x 2 )<p s K ,(x 3 ,b 3 )}E f (t 2 a )h a (x 3 ,x 2 ,b 3 ,b 2 )], (A5) 

\ rl roo 

M 3 = —==<6AtxC f M% / dxidx 2 dx 3 \b\db\b 2 db 2 <j>B(xi,b 1 )(f) D (x 2 ,b 2 ) 
y2JS c Jo Jo 

x \{(x 3 + (x 2 - 2x 3 - l)rj)(j) K *(x 3 ,b 2 ) + r 2 r 3 (2 + x 2 + x 3 ) (f) s K *(x 3 , h) 

-r 2 r 3 (x 2 - x 3 ) (Pk,(x 3 , b 2 )}E m (t 1 m )h^\x 1 ,x 2 , x 3 ,h, b 2 ) 
-{x 2 (p K *(x 3 ,b 2 ) + r 2 r 3 (x 2 + x 3 ) (p s K ,(x 3 , b 2 ) 

+r 2 r 3 (x 2 - x 3 )(j) t Kf (x 3 ,b 2 )}E m (t 2 m )h (2) (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) . (A6) 
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d. B+ -► D* + K*° decay 



F L = lQirC F M\ 



rl roo 

j dx 2 dx 3 / b 2 db 2 b 3 db 3 4>d*(x 2 , b 2 ) 
Jo Jo 



x [{(x 3 - 2x 3 r^ + r^)0 x *(a;3,63) + r 2 r 3 (^*(:r 3 , 63) - ^^.(xs, 6 3 )}S/(t„)/i (x2,a:3, 62, 63) 



-{(1 - r^)x20^.(x 3 ,6 3 ) + 2r 2 r 3 (x 2 - l)<f) s K *(x 3 ,b 3 )}E f (t 2 a )h a (x 3 ,x 2 ,b 3 ,b 2 ) 



(AT) 



M L = —==6471(7^^^1 dx!dx 2 dx 3 \ b\db\b 2 db 2 <f)B(xi,b 1 )(f) D *(x 2 ,b 2 ) 
y2JS c Jo Jo 

x [{(-x 3 + (x 2 + 2x 3 - l)?^)^*^,^) + r 2 r 3 (x 2 - x 3 ) <p s K *(x 3 ,b 2 ) 
+r 2 r 3 (2-X2- x 3 ) 0^*(x 3 , b 2 )}E m (t 1 m )h£\x 1 ,x 2 , x 3 , 61, 62) 
+{ (l - 2rf ) x 2 0^* (x 3 , 6 2 ) + r 2 r 3 (^2 - x 3 ) <f) s K * (x 3 , b 2 ) 

+r 2 r 3 (x 2 + x 3 ) 0x-.(x 3 , 6 2 )}£'m(0 /i i 2) ( a; i ) ^2, £3, 61, 62) • 



(A8) 



Fti — 



M' 



Tl 



rl roo 

l6irC F Mg dx 2 dx 3 b 2 db 2 b 3 db 3 (f)D*(x 2 ,b 2 ) 
Jo Jo 

x [{-2r 2 r 3 x 3 (j) v K * (x 3 , b 3 ) - 2r 2 r 3 x 3 (f> a K * (x 3 , 63) + 2r\f K * (x 3 , b 3 )}E f (tl)h a (x 2 , x 3 , b 2 , b 3 ) 
+{2r 2 r 3 0^ (x 3 , b 3 ) + 2r 2 r 3 <j) a K * (x 3 , b 3 )}E f (t 2 a )h a (x 3 , x 2 , 63, b 2 )] , (A9) 

: /7T1 r^ 647rCFM^ / dxidx 2 dx 3 / b x db x b 2 db 2 (j)B(x 1 ,b 1 )(f) D *(x 2: b 2 ) 
\/2JS c Jo Jo 

: {2r 2 r 3 <f) v Kt (x 3 , b 2 ) + 2r 2 r 3 ^. (a*, b 2 ) - 2r\ (1 - x 2 ) f K * (x 3 , b 2 )} 
E m (t 1 m )h^ ) (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) - {2rjx 2 (i)K t (x 3 ,b 2 )}E rn (t 2 m )h ( a\x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) . (A10) 



rl roo 

F T2 = VoixC F M\ \ dx 2 dx 3 \ b 2 db 2 b 3 db 3 (j)D*(x 2 ,b 2 ) 
Jo Jo 

x {-2r 2 r 3 0^(x 3 ,6 3 ) +2r 2 r 3 (f) a K *(x 3 ,b 3 )}E f (t 1 a )h a (x 2 ,x 3 ,b 2 ,b 3 ) 
+{2r 2 r 3 x 2 (p v K *(x 3 , b 3 ) - 2r 2 r 3 x 2 0x*( a; 3, h)}E f (t 2 a )h a (x 3 , x 2 , b 3 , b 2 ) 



(All) 



M T 2 = 



<oAixC f Mb I dx\dx 2 dx 3 \ b\db\b 2 db 2 (f>B(xi,bi)<f) D *(x 2 ,b 2 ) 
Jo Jo 



V2N, 



x {2r 2 r 3 <p v K *(x 3 ,b 2 ) - 2r 2 r 3 (j) a Kt (x 3 ,b 2 )}E m (t 1 m )h£\x 1 ,x 2 ,x 3 ,b l ,b 2 ). 



(A12) 



where the subscript L(T(1,2)) , i.e. the helicity states of the two vector mesons o in 
eq. (17), stands for the longitudinal (transverse) component respectively. Conveniently 
we choose the polarization state Tl as e 2F = ^(0, 0, 1, i), e 3 T = ^75 (0, 0, 1, — i), T2 as 
£2T = ^(0, 0, 1, — i),e 3 r = ^7= (0,0, l,i). Each amplitude A a also is the sum of two parts, 
factorizable and non-factorizable diagrams, related by eq. (16). 
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2. The ifM meson wave functions 



The K and K* meson wave functions are given as [19,20] 

i a i __\ Ik 



4> A K {x) 



2^2N ( 

f. 



Sx(l - x) {l + 0.51(1 - 2x) + 0.2C 2 3/2 (1 - 2x)} , 



<Pk(x) = ^== {l + 0.212C 2 1/2 (1 - 2x) - 0.148C 4 1/2 (1 - 2x)} 



t>K*(X 



-^=Qx(l -x)\l + 0.57(1 - 2x) + 0.07C 2 3/2 (1 - 2x) 
y2N c L 



K 



2 v / 2iVc 



1 - 2x) {l + 0.1581[-3 + 5(1 - 2a;) 2 ]} , 



(t) l K ,{x) = 2 ^ r {°- 3 ( 1 - 2a; ) [ 3 ( X - 2a; ) 2 + - 2x ) 1-68C 4 1/2 (1 - 2x) 

+0.06(1 - 2xf [5(1 - 2xf - 3] + 0.36 [1 - 2(1 - 2x)(l + ln(l - x))\ 
f. 



<t>K*{x) 



3(1 -2x) 1 + 0.2(1 - 2a:) + 0.6(10a; 2 - 10a: + 1) 



2V2N C 

-0.12x(l -x) + 0.36[1 - 6a: - 2 ln(l - , 



2v/2iV; 



6x(l - x) 1 + 0.6(1 - 2x) + 0.04C 2 3/2 (1 - 2x) 



MO I \ f K * f 3 



[X = 



1 + (1 - 2xf + 0.44(1 - 2xf\ + 0.4C 2 1/2 (1 - 2x) 
+0.88C 4 1/2 (1 - 2x) + 0.48[2a; + ln(l - , 

f ' ' {3(1 - 2a;) [l + 0.19(1 - 2x) + 0.81(10:r 2 - 10a: + 1) 



-1.14x(l -x) + 0.48[1 - 6a: - 2 ln(l - x)]} , 
with the Gegenbauer polynomials, 

Cl'\0 = ^ - 1) , Q 1/2 (O = ^(35£ 4 -30£ 2 + 3), Q 3/2 (0 = ^(5e 2 -l). 



(A13) 
(A14) 
(A15) 

(A16) 
(A17) 

(A18) 
(A19) 

(A20) 

(A21) 
(A22) 



3. Some used formulae 

The definitions of some functions used in the text are presented in this appendix. In the 
numerical analysis we use one loop expression for strong coupling constant, 



a s (n) 



An 



(A23) 
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where (3q = (33 — 2n/)/3 and rif is number of active flavor at appropriate scale. A is QCD 
scale, which we use as 250 MeV at nj = 4. We also use leading logarithms expressions for 
Wilson coefficients Ci )2 presented in ref. [15]. 

The function Ef and E m including Wilson coefficients are defined as 



E f (t) =a(t)a s (t) e -s D (t)-s K (t)^ 
E m (t) = Ci(*K(t) e s B (t)-s D (tys K (t) 



b^=b2 , 



where 



a(t) = C 2 (t) + 



Ci(t) 



(A24) 
(A25) 



(A26) 



and Sb, S b , and Sk result from summing both double logarithms caused by soft gluon 
corrections and single ones due to the renormalization of ultra-violet divergence. The above 
Sb d k are defined as 



S B (t) = s(x 1 P+,b 1 )+2 f % lq {»'), 

ft d u > 

S D (t) = s(x 2 P 2 + ,b 2 ) + 2 / -%(/*'), 

J 1/62 fJ> 



S K (t) = s(x 3 P 3 ,b 3 ) + s((l-x 3 )P 3 ,b 3 ) + 2 f % lq {n'), 



where s(Q,b), so-called Sudakov factor, is given as [24] 



1/63 fi 



s( Q, 6 ) = /^[{2 (27E _ 1 _ ln2) + CflD ^Wrt 
Jl/bfJ,' [{3 //J 71 



(A27) 
(A28) 
(A29) 



+ 



67_ TT^ 
9 3 



IE 



10 2 , e 
+ In — 



'a s (jj!) 



7T 



(A30) 



7£ = 0.57722 • • • is Euler constant, and 7 9 = —a s /ir is the quark anomalous dimension. 

The functions h a , and in the decay amplitudes consist of two parts: one is the 
jet function S t (xi) derived by the threshold resummation [14], the other is the propagator 
of virtual quark and gluon. They are defined by 

K(x 2 ,x 3 ,b 2 ,h) = S t (l-x 3 ) [^-j H^\M B ^(l-rl)x 2 x 3 b 2 ) 

x \h£\m b ^(1 - r 2 2 )x 3 b 2 )J (M B ^(l-r 2 2 )x 3 b 3 )6(b 2 - b 3 ) + (6 2 ^ 63)} , (A31) 

(A32) 



h { J ) (x 1 ,x 2 , x 3 , 61, 6 2 ) = {f 4 1) (M BV /(1 -ri)x 2 x 3 b 1 )Jo(M BV /(l -ri)x 2 x 3 b 2 )^(b 1 - b 2 ) + (bi <-> b 2 ) 



x 



KoiMsF^h), for 



2 



^(MbJIFXJ bi), for F ( 2 ?) <0 



F U) > 



(j) 



(?) 
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where Hq^(z) = Jo(z) + i Y (z), and F^s are defined by 

=x 1 + x 2 + (l-x 1 - x 2 )x 3 (l - r 2 2 ), (A34) 
F^=x 3 ( Xl -x 2 )(l-r 2 2 ). (A35) 

We adopt the parametrization for S t (x) of the factorizable contributions, 

S t (x)^ ^^^±^- [x(l-x)] c , c = 0.3. (A36) 

In the non-factorizable annihilation contributions, St(x) gives a very small numerical effect 
to the amplitude [14]. Therefore, we drop S t (x) in and The hard scale t's in 

the amplitudes are taken as the largest energy scale in the H to kill the large logarithmic 



radiative corrections: 

t\ = max(M BV /(l-ri)x 3 , l/b 2 , l/b 3 ), (A37) 

t 2 a = max(M B ^(1 - r|)x 2 , l/b 2 , l/b 3 ), (A38) 

iP m = max(M Bv ^Fy , - r|)x 2 x 3 , 1/bi, l/b 2 ). (A39) 
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